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We show that it is possible to perform a continuous measurement that continually projects a 
nano-resonator into its energy eigenstates by employing a linear coupling with a two-state system. 
This technique makes it possible to perform a measurement that exposes the quantum nature of the 
resonator by coupling it to a Cooper-pair Box and a superconducting transmission-line resonator. 
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It is now possible to construct nanomechanical res- 
onators with frequencies on the order of 100 MHz, and 
quality factors of 10 5 [1-9]. This opens up the excit- 
ing prospect of observing quantum behavior in meso- 
scopic mechanical systems, implementing quantum feed- 
back control in these devices [10, 11], and exploiting 
them in technologies for such applications as metrology 
and information processing [12]. The position of these 
resonators can be monitored by using a single electron 
transistor (SET) placed nearby [13-15], and such a mea- 
surement has recently been realized close to the quan- 
tum limit by Schwab et al. [7, 9]. However, to observe 
the quantum nature of a nano-resonator one must mea- 
sure an observable that is not linear in the position or 
momentum of the resonator, and such measurements are 
considerably more difficult to devise. One approach that 
has been investigated is to perform a Quantum Non- 
Demolition (QND) measurement of the resonator energy 
which would project the resonator into its (discrete) en- 
ergy eigenstates. This would result in the observation of 
jumps between these states, a clear signature of quantum 
behavior. However, such a measurement requires the con- 
struction of a non-linear interaction with the resonator, 
and devising such a coupling with sufficient strength is 
challenging [16, 17]. 

Here we show that it is possible to perform a mea- 
surement that continually projects the resonator onto the 
basis of energy eigenstates (the Fock basis) using only a 
linear coupling to a probe system. While the resulting 
measurement is not a QND measurement, it neverthe- 
less allows a direct observation of the quantum nature of 
the resonator because it continually projects the system 
onto the Fock basis. The method exploits the fact that 
the linear coupling will transfer the effect of a non-linear 
measurement of the probe onto the resonator, and has 
similarities with that used in atom optics in the detection 
of quantum jumps with resonance fluorescence [18-21]. 
The measurement technique we describe is also expected 
to have applications to state-preparation and feedback 
control [11]. 



Before we begin we briefly discuss the anatomy of a 
quantum measurement. To perform a measurement of 
an observable A of a quantum system one couples the 
system to a second "probe" system. If one choses an 
interaction Hamiltonian H = hXAB, where B is an ob- 
servable of the probe, then after a time t, this will cause 
a shift of XAt in the probe observable conjugate to B, 
which we will call C. This shift in C can be measured 
to obtain the value of A. The observable B is chosen so 
that its conjugate observable C can be easily measured 
directly by an interaction with a classical apparatus. To 
obtain a continuous measurement of A one proceeds in 
an analogous fashion, except that the interaction is kept 
on, and C is continually monitored. Such a measurement 
provides a continual stream of information about A, and 
is usually referred to as a continuous measurement [22]. 
Such a measurement will continually project the system 
onto the eigenstates of A. The measurement is referred to 
as a QND measurement if A commutes with the system 
Hamiltonian, so that the system remains in an eigenstatc 
of A once placed there by the measurment [23] . 

We now consider coupling a nano-resonator to a sec- 
ond harmonic oscillator via an interaction linear in the 
resonators position: H = HXxpXp. Here xp = a + a' 
is the resonator position, and xp = b + is the posi- 
tion of the probe oscillator which we will take to have 
the same frequency as the resonator, wr. This coupling 
transfers energy between the two oscillators, as well as 
correlating their phases. If A <C then the rotating 
wave approximation gives H = hX(ab^ + baJ) which is an 
explicit interchange of phonons. 

Now consider what happens if we perform a continu- 
ous measurement of the energy of the probe oscillator. 
(This would be a QND measurement of the probe if it 
were not coupled to the resonator.) Since the energy of 
the resonator continually feeds into the probe oscillator 
(and vice versa), this measurement of the probe must tell 
us about the energy of the resonator. We should there- 
fore expect the measurement to localize both the probe 
and the resonator to their energy eigenstates. This is 
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somewhat surprising from the point of view of the dis- 
cussion above, since the (linear) position coupling would 
be expected to transfer phase information to the probe, 
disturbing the energy eigenstates and projecting instead 
onto the position basis. The two arguments can be rec- 
onciled by noting that the phase information regarding 
the system is contained in the phase of the probe, and 
this phase is continually destroyed by the energy mea- 
surement on the probe. Since the probe is projected into 
an energy eigenstate, the interaction does not imprint a 
phase back on the resonator from the probe, and there is 
nothing to prevent it localizing to an energy eigenstate. 
Nevertheless, the expected disturbance to the Fock states 
is not completely eliminated, as we shall see, because the 
linear interaction causes jumps between these states. 

To test the above intuition, we now simulate the evo- 
lution of the coupled oscillators. The stochastic master 
equation (SME) describing their dynamics is 

dp = -(i/h)[H,p]dt-k[N,[N,p]}dt 

+Ak[Np + pN -2(N)p](dr- (N)dt), (1) 

where N = b^b is the phonon number operator for the 
probe, H = K\{aw +baJ), and k is the strength of the en- 
ergy measurement on the probe. The observers measure- 
ment record is dr = (N)dt+dW/ VSk, where dW is Gaus- 
sian white noise satisfying the relation dW 2 = dt [24]. 
The observer obtains p(t) by using her measurement 
record to integrate Eq.(l). The simulation is performed 
using a second order integrator for the deterministic mo- 
tion, and a simple half-order Newton integrator for the 
noise term. This involves picking a random Gaussian 
variable with variance At at each time-step At. We 
choose the initial states of the two oscillators as coher- 
ent states with mean phonon number 2, and measure 
time in units of k. Since there is no additional noise ap- 
part from that induced by the measurement, we can use 
the stochastic schrodingcr equation equivalent to Eq.(l), 
which reduces the numerical overhead [25]. 

The results of the simulation are depicted in Figure 1. 
Setting the initial interaction strength to A = fc/20, we 
find as expected that the resonator's energy variance re- 
duces essentially to zero at rate A. The measurement 
thus projects the system onto the energy eigenbasis as 
required. If we start the probe system in a known 
energy eigenstate (by measuring its energy before we 
switch on the interaction), then the measurement pro- 
cess also provides full information regarding the initial 
energy of the resonator, as required of a measurement 
of energy. However, the interaction causes an additional 
effect: the two oscillators undergo equal and opposite 
quantum jumps between their energy eigenstates. (Af- 
ter a time of t = 50/fc we reduce A. This reduces the 
rate of jumps so that both the jumps and the periods of 
stability are clearly visible.) This behavior can be under- 
stood as follows. The energy measurement tends to keep 
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FIG. 1: Here we plot the evolution of the nanomechanical 
resonator and the probe oscillator: (a) the energy of the res- 
onator; (b) the energy of the probe; (c) the variance of the 
energy of the resonator. The interaction strength A = fc/20 
for tk < 50, and A = 7.5 x 10~ 3 fc for tk > 50. 

the resonator and the probe in their energy eigenstates 
because of the quantum Zeno effect. However, the inter- 
action is continually trying to transfer energy between 
the two oscillators, and at random intervals this over- 
comes the quantum Zeno effect and the two oscillators 
jump simultaneously between energy states. The jumps 
are equal and opposite and thus preserve their combined 
phonon number. The rate of the jumps is determined 
by the relative size of A and k: when k 3> A the jumps 
are suppressed by the quantum Zeno effect, the energy 
transfer rate is reduced, and correspondingly the rate of 
information extraction from the system is also reduced. 

So far we have been considering a Harmonic oscillator 
as the probe system. We note now that since a harmonic 
oscillator truncated to it lowest two energy levels is a 
two-level system, this suggests that one might be able to 
use a two-level system as a probe in the same way. This 
would increase considerably the range of possible experi- 
mental realizations. We find that this is indeed the case; 
a two-level system (TLS) is similarly effective at project- 
ing the nano-resonator onto an energy eigenstate. If we 
truncate the probe harmonic oscillator to its lowest two 
levels, then the interaction between the system and probe 
becomes H = Act^xr, and the energy measurement on 
the probe is a measurement of a z . In figure 2(a), we plot 
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FIG. 2: Here we plot the evolution of the nanomechanical 
resonator under a measurement by a Cooper-pair box, as well 
as histograms of the distribution of the average phonon num- 
ber over the evolution, (a) Zero temperature (T = 0) and 
no damping (r = 0); (b) T = 6mK and T = fc/500; (c) 
T = 32mK and F = fc/2500. 



the evolution of the energy of the resonator under such 
a measurement. Whereas in our previous simulations we 
assumed that A -C wr and made the rotating wave ap- 
proximation, here we make no approximation. We choose 
k = ur/20, and A = (3/4)fc, so that there are rapid ex- 
changes of energy between the two systems. Once again 
the variance of the nano-resonator's energy reduces as 
required, but this time the resonator jumps (rapidly) be- 
tween only two adjacent energy levels, since the TLS has 
only two energy states. We also see a new effect due to 
the fact that the total number of excitations is no longer 
preserved by the interaction (because we have not made 
the rotating wave approximation). Because of this the 
resonator gets energy kicks from the TLS that are not 
associated with a flip of the TLS state. These are occa- 
sionally sufficient to cause the resonator to jump between 
phonon states, shifting the offset of the rapid oscillations 
up or down by one phonon. 

We now show how such an energy measurement can 
be implemented using a Cooper-Pair-Box (CPB) cou- 
pled in turn to a superconducting transmission-line res- 
onator [26]. A CPB is a superconducting island, whose 
two charge states consist of the presence or absence of a 
Coopcr-pair on this island. If we work at the degeneracy 
point where the two charge states have the same charging 



energy, then the free Hamiltonian of the CPB contains 
only the Joscphson tunneling term Eja x . If we place 
a bias voltage with frequency 8 on the nano-resonator, 
and place the CPB adjacent to it, we obtain the coupling 
term Acos(<5i)cr z XR [27]. If we then place the CPB in a 
superconducting resonator (SR), and detune the Joseph- 
son tunneling frequency Ej / h from the SR frequency cos 
by an amount A, then the interaction between the CPB 
and the superconducting resonator is well approximated 
by the Hamiltonian H = h(g 2 / A)<r x c < c [28]. Here c is 
the annihilation operator for the SR mode and g is the 
so-called "circuit QED" coupling constant between the 
CPB and the SR [29] . This approximation requires that 
{Ej/h, ujs) ^> A and A ^> g, and we set 8 = Ej /H—ljr to 
bring the resonator-CPB interaction on-resonance. Thus 
the full Hamiltonian for the nano-resonator, the CPB and 
the superconducting resonator is 
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This achieves the required configuration in which the 
nano-resonator is coupled to a CPB via one Pauli op- 
erator, and the CPB is coupled to a second probe system 
via a second Pauli operator. All that is required now is 
that we use the second probe system (the SR) to per- 
form a continuous measurement of a x . The interaction 
term h(g 2 / A)a x c^ c means that the a x eigenstates of the 
CPB generate a frequency shift of the SR, which in turn 
produces a phase shift in the electrical signal carried by 
a transmission line connected to the SR. Two methods 
for continuously monitoring this phase shift with high fi- 
delity have been devised by Sarovar et al. [28]. If one 
performs a continuous measurement of the phase of the 
SR output signal, then one can adiabatically eliminate 
the SR and obtain an equation describing the continuous 
measurement of the CPB (this type of adiabatic elimina- 
tion procedure is detailed in [30]). The resulting SME is 
precisely Eq.(l), where the Hamiltonian H is replaced by 
H and the phonon number N is replaced by a x . The im- 
portant quantity is the final measurement strength k of 
this a x measurement. The adiabatic elimination results 
in the measurement strength k = (g 4 \a\ 2 )/(A 2 j). where 
7 is the decay rate of the SR, and |a| 2 is the average num- 
ber of photons in the SR during the measurement. The 
adiabatic elimination requires that 7 ^> g 2 /A. We note 
that this second inequality is merely required to ensure 
the accuracy of our expression for k — the measurement 
can be expected to remain effective without it. 

Readily obtainable values for the circuit QED pa- 
rameters are g = it x 10 8 s _1 and 7 as low as 6 x 
10 6 s _1 [29]. A realistic frequency for the nano-resonator 
is ujr/2it = 100 MHz and for the superconducting res- 
onator is ujs/Z-k = 10 GHz. We choose the CPB fre- 
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4tt x 10 s s _1 and set g = A/40, 
50. With 7 = 7T x 10 7 s" 1 we 
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then have 'jA/g 2 = 20. With these parameters, choosing 
even a modest value of \a\ 2 — 2 x 10 3 provides a mea- 
surement strength of k = 4 x 10's _1 . The interaction 
strength A is not a limitation, and can easily be as high as 
10 8 s _1 [11, 27]. Nano-resonators typically have quality 
factors of Q = 10 5 , giving a damping rate of T w 10 4 s _1 . 

We now turn to the question of observing the quantum 
nature of the resonator. In this measurement scheme 
the quantum behavior of the resonator is not manifest 
in energy jumps resulting from exchanges of excitation 
number with the CPB; even if the resonator were clas- 
sical these jumps would occur because the CPB states 
are discrete. The discrete nature of the resonator states 
are manifest in the fact that the energy measurement lo- 
calizes the resonator energy to integer multiples of Swr, 
rather than just any value consistent with the thermal 
distribution. As a result the rapid oscillations due to 
excitation exchanges only occur between these discrete 
values (to within the energy localization induced by the 
measurement). Further, thermal noise does not cause 
the oscillator to undergo Brownian motion as it would 
during a continuous energy measurement on a classical 
oscillator, but instead induces quantum jumps between 
the discrete levels. As a result a histogram of (N) over 
time is therefore peaked at integer values, in sharp con- 
trast to the classical case. 

Since we can achieve k 3> T, we would expect to be 
able to observe the discreteness of the energy levels at 
low temperatures. We now perform numerical simula- 
tions to verify this. These simulations are numerically 
intensive, so we make the rotating wave approximation, 
and to include the thermal noise we use an approxima- 
tion to the Brownian motion master equation [31] that 
takes the Lindblad form [32]: p = r(£ + l)D[o/2]p + 
r£X>[aT /2]p. Here T is the damping rate of the resonator, 
£ = coth[fkon / (JlksT)] (where T is temperature), and 
T>[c\p = 2cpS —Sep— pS c for any operator c. The CPB 
is also subject to damping and dephasing, and we include 
both of these at a rate 7cpb = 1 x 10 6 s _1 , which is not far 
from current values [33]. We plot the results in Figure 2 
using the parameters given above with A = (3/4)/:. Fig- 
ure 2(b) shows the results for T = fc/500 and T = 6 mK 
and Figure 2(c) for T = fc/2500 and T = 32 mK. For each 
case we plot the histogram of (TV), and this shows that 
the peaks at integer values are clearly visible. We also 
see that the effect of the thermal noise is larger when the 
resonator is in higher energy eigenstates; as the phonon 
number increases the peaks are washed out and the be- 
havior of (N) becomes indistinguishable from Brownian 
motion. 
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